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THE  CONSISTENT  SECOND-ORDER  THEORY 
OF  WAVE/STRUCTURE  INTERACTION 

by 
C.  J.  Garrison 


ABSTRACT 

The  consistent  second-order  theory  of  the  interaction  of  regular  gravity 
waves  with  a  fixed  object  in  water  of  finite  depth  is  developed.  The  theory  is 
carried  out  for  the  most  general  case  of  a  body  of  arbitrary  shape  which  may 
extend  through  the  free-surface  or  be  completely  immersed.  The  incident  wave 
evolves  in  the  development  as  a  second-order  Stokes'  wave.  Boundary-value 
problems  are  established  for  both  the  first-  and  second-order  velocity  poten- 
tials and  a  numerical  method  based  on  the  Green's  function  is  outlined. 

INTRODUCTION 

The  determination  of  forces  exerted  by  gravity  waves  on  large  structures 
immersed  in  the  sea  has  become  of  great  practical  interest  in  recent  years.  For 
example,  in  the  design  of  bottom-mounted  oil  storage  facilities  or  large  ocean 
caissons,  the  wave-induced  horizontal  and  up-lift  forces  and  overturning  moments 
are  factors  of  primary  importance.  The  effect  of  large  amplitude  waves  in  parti- 
cular is  of  importance  in  the  determination  of  the  permanence  of  an  ocean  struc- 
ture and,  therefore,  a  higher-order  theory  appears  to  have  significant  practical 
value.  For  example,  Apelt  and  Macknight  (1976)  found  measured  forces  on  a  ocean 
caisson  model  in  fairly  large-amplitude  shallow-water  waves  to  be  considerably 
in  excess  of  calculations  based  on  linear  diffraction  theory. 

With  this  application  in  mind  then,  the  solution  through  the  second-order 
in  wave  height  is  developed  herein  for  regular  wave  interaction  with  a  fixed 
object  of  arbitrary  shape  in  water  of  finite  depth.  The  theoretical  development 


is  carried  out  in  a  mathematically  consistent  manner  in  that  all  terms  through 
the  second-order  in  wave  height  are  included  in  the  analysis,  and,  as  such, 
represents  the  extension  of  the  linear  theory  presented  previously  by  Garrison 
and  Seetharama  Rao  (1971).  The  problem  is  recognized  as  a  regular  perturbation 
problem  in  the  small  parameter  (wave  height/characteristic  body  dimension)  and 
the  incident  wave  appears  within  this  framework  as  the  second-order  Stokes' 
wave.  A  numerical  method  for  solving  the  resulting  equations  suitable  for  digital 
computer  evaluation  is  outlined. 

REVIEW  OF  LITERATURE 

The  solution  of  the  linear  wave/structure  interaction  problem  is  now  fairly 
well -developed  for  infinite  depth  as  well  as  for  the  finite  depth  case.  A 
number  of  papers  have  appeared  in  the  literature  dealing  with  the  oscillation  of 
two-dimensional  bodies  on  a  free  surface  in  water  of  infinite  depth,  a  problem 
which  is  mathematically  similar  to  the  fixed-body/wave  interaction  problem. 
Examples  include  the  work  of  Ursell  (1949),  Porter  (1960),  Vugts  (1968),  and 
Paul  1 ing  and  Richardson  (1962).  Dean  and  Ursell  (1959)  made  both  an  experimental 
and  theoretical  study  of  small  amplitude  wave  interaction  with  a  fixed,  semi- 
submerged  circular  cylinder  in  deep  water.  Yu  and  Ursell  (1961)  treated  the 
problem  of  a  semi -submerged  circular  cylinder  oscillating  vertically  in  water  of 
finite  depth. 

Somewhat  less  attention  has  been  given  to  solving  the  corresponding  three- 
dimensional  problems.  Apparently  the  first  analysis  of  the  linear  wave  inter- 
action with  a  fixed  body  was  carried  out  by  Havelock(1940)  for  the  case  of  a  pile 
in  water  of  infinite  depth  and  MacCamy  and  Fuchs  (1954)  solved  the  same  problem 
for  water  of  finite  depth.  Havelock  (1959)  also  evaluated  the  added  mass  and 
damping  coefficients  for  a  heaving  semi -immersed  sphere  in  deep  water,  and  Wang 
(1966)  extended  this  to  include  the  finite  depth  case.  Haskind's  relations 
(Haskinds  (1957))  as  discussed  by  Newman  (1962)  may  be  used  to  determine  the 


linear  or  first-order  forces  acting  on  the  fixed  body  from  a  knowledge  of  the 
damping  coefficients  for  the  same  body  oscillating  in  still  water.  Garrison  (1974) 
has  made  this  evaluation  using  Havelock's  results  for  the  hemisphere  and  the  re- 
sults were  in  good  agreement  with  his  direct  calculation  of  the  heave  force  based 
on  the  Green's  function  approach.  Kim  (1965,  1966)  has  applied  the  Green's 
function  approach  to  compute  the  excitation  forces,  added  mass  and  damping  co- 
efficients for  a  semi -immersed  ellipsoid  in  deep  water. 

Most  work  on  wave/body  hydrodynamics  reported  in  the  literature  treats  con- 
figurations and  conditions  which  have  primary  application  to  the  study  of  ship 
motion;  only  recently  an  interest  in  large,  fixed  offshore  structures  has  develop- 
ed. Even  though,  probably  the  first  potential  flow  (or  diffraction)  solution  to 
a  wave  force  problem  having  direct  application  to  large  ocean  structures  was  ob- 
tained over  twenty  years  ago  by  MacCamy  and  Fuchs  (1954),  namely,  the  small 
amplitude  wave  interaction  with  a  fixed  vertical  circular  cylinder  (pile)  in 
water  of  finite  depth.  This  diffraction  solution  represents  the  only  closed  form 
solution  of  its  type  available  in  the  literature  and,  therefore,  is  of  primary 
importance  to  this  general  problem.  Little  work  of  this  type  having  application 
to  large  ocean  structures  has  appeared  since  that  time  until  recently  when 
Garrison  and  Seetharama  Rao  (1971)  applied  the  Green's  function  approach  to  calcul- 
ate the  first-order  pressure  distribution  and  resulting  forces  acting  on  a  bottom- 
mounted  hemisphere.  Their  experimental  results  for  this  configuration  compared 
well  with  the  theory.  Garrison  and  Chow  (1972)  then  extended  this  analysis  to 
include  fixed  bodies  of  arbitrary  shape  and  compared  their  theoretical  results 
with  experimental  results  corresponding  to  two  different  submerged  oil  storage 
vessel  configurations.  Milgram  and  Halkyard  (1971)  also  have  developed  a  linear 
theory  and  applied  their  method  to  certain  axisymmetric  bodies  in  deep  water. 
Garrison  (1974,  1977)  has  more  recently  given  details  of  a  practical  method 
based  on  the  Green's  function  to  evaluate  added  mass  and  damping  coefficients  for 
floating  bodies  of  arbitrary  shape  in  water  of  finite  depth,  and  at  the  same  time 


calculate  the  linear  wave  forces. 

It  may  be  stated  that,  in  general,  linear  theory  is  quite  well  developed 
and  good  agreement  between  theory  and  experiment  has  been  obtained  for  wave  forces 
corresponding  to  small  amplitude  waves  in  both  the  two-  and  three-dimensional 
cases.  This  is  true  as  well  for  bodies  oscillating  in  a  free  surface  in  still 
water.  Experience  has  also  shown  that  linear  theory  gives  good  wave  force  results 
when  the  body  is  fairly  deeply  submerged  even  if  the  amplitude  of  the  incident 
wave  is  not  small.  However,  if  the  object  is  surface  piercing  or  not  deeply  sub- 
merged and  the  water  is  shallow,  nonlinear  effects  caused  by  waves  of  finite  amp- 
litude become  pronounced.  In  such  cases  linear  theory  is  inadequate,  and  a  higher 
order  wave/structure  interaction  theory  becomes  necessary. 

Nonlinear  potential  solutions  to  wave/structure  interaction  problems  are 
much  more  limited  than  linear  solutions.  Olgilvie  (1963)  solved  for  the  first- 
and  second-order  forces  on  a  horizontal,  submerged  circular  cylinder  in  water  of 
infinite  depth,  the  wave  crests  being  parallel  to  the  cylinder  axis.  However, 
his  second-order  forces  included  only  the  time-average  part  and,  as  a  consequence, 
the  first-order  potential  only  was  needed  and  obtained.  This  work,  therefore,  did 
not  actually  represent  a  complete  second-order  solution.  However,  a  consistent 
second-order  theory  for  a  horizontal  cylinder  in  still  water  has  been  developed 
by  Lee(1968)  and  Potash(1970) .  Garrison  and  Smith(1977)  have  treated  the  corres- 
ponding wave/fixed-body  interaction  problem. 

Recently  there  has  been  some  interest  in  developing  analyses  for  nonlinear 
wave  interaction  with  three  dimensional  fixed  bodies  in  water  of  finite  depth. 
Chakrabarti  (1972)  made  an  attempt  at  the  extension  of  MacCamy  and  Fuch's  (1954) 
linear  solution  for  the  vertical  circular  cylinder  (pile)  to  the  fifth  order. 
However,  he  treated  the  free-surface  boundary  condition  improperly  so  that  the 
results  are  mathematically  inconsistent  and,  consequently,  of  little  value  for 
making  judgements  regarding  the  magnitude  of  nonlinear  diffraction  effects.  More 
recently  Raman,  et.al.  (1976,  1977)  have  also  treated  this  problem 


although  the  method  of  solution  is  extremely  complex  and  there  appears  to  be  no 
simple  means  of  checking  the  results. 

Isaacson  (1977)  has  argued  that  the  second-order  solution  for  a  body  with 
a  circular  waterline  curve  does  not  exist  because  the  free-surface  and  kinematic 
boundary  condition  on  the  cylinder  are  not  consistent  at  the  mean  waterline  on 
the  cylinder.  However,  this  inconsistency  does  not  preclude  the  existance  of  a 
solution;  it  appears  that  such  inconsistencies  at  the  point  where  two  different 
boundary  conditions  join  is  rather  common. 

In  this  report  the  solution  through  the  second-order  in  wave  height  is 
developed  for  bodies  of  arbitrary  shape  held  fixed  in  regular  waves  in  water  of 
finite  depth.  The  problem  is  formulated  as  a  regular  perturbation  problem,  and 
the  incident  wave  is  shown  to  represent  a  second-order  Stoke1 s  wave.  It  is  also 
shown  that  the  second-order  problem  is  the  same  type  of  boundary-value  problem 
as  the  first-order  problem,  the  primary  difference  being  that  the  second-order 
free  surface  boundary  condition  is  nonhomogeneous.  The  solution  is  formulated 
for  both  the  first-  and  second-order  potentials  in  terms  of  a  Green's  function, 
and  the  numerical  results  are  presented  for  several  configurations. 


FORMULATION  OF  THE  PROBLEM 

The  problem  under  consideration  is  depicted  in  Figure  1.  A  rigid  object 
having  a  characteristic  dimension  a  is  immersed  in  water  of  depth  h,  and 
a  train  of  regular  waves  propagates  in  the  positive  x-direction.  The  poten- 
tial solution  is  sought  to  the  interaction  of  the  incident  wave  with 


^(x,z,  t) 


S(x,y,  z  )=  0 
(Surface  of  body) 


y=  0   plane :   mean  water  elevation 


FIGURE  1   DEFINITIONS 


the  fixed,   rigid  body  through  the  second-order  in  the  wave  amplitude, 
or  more  precisely,   in  the  ratio  of  the  wave  amplitude  to  characteristic 
dimension  of  the  object,  a". 

The  fluid  is  assumed  to  be  irrotational   so  that  the  velocity 
potential,   *,  may  be  defined  as 

where  q~  denotes  the  fluid  velocity  vector.     Assuming  the  fluid  to  be 
incompressible,  it  follows  that  the  velocity  potential  must  satisfy 

VZ$Cx/y,z/£)  =  o  (2) 

within  the  fluid  region.  (The  barred  quantities  denote  dimensional 
quantities. ) 

On  the  bottom  the  kinematic  condition  specifying  zero  normal 
velocity  is 

$q(*,-h,2,t)   =0  (3) 

where  h  denotes  the  mean  fluid  depth.     On  the  rigid,  wetted  surface  of 
the  body,  specified  by      SCx,Q,l)  -    O  »   the  kinematic  boundary 

condition  describing  the  zero  normal   velocity  condition  is  given  by 

Y<§>  '    V5   =o  (4) 

The  elevation  of  the  free  surface  above  the  y  =  0  plane  is  denoted 
by  5"(x,z,t)  and  the  kinematic  boundary  condition, 


Pt 


=  o 


(5) 


must  be  satisfied  on  y  =  n.  Moreover,  on  the  free  surface  the  pressure 
is  set  equal  to  zero,  and,  accordingly,  Bernoulli's  equation  supplies  the 
dynamic  free  surface  boundary  condition, 

^.c*,?,^?;^/^^  (6) 

where  h  denotes  the  Bernoulli  constant. 

It  is  convenient  to  cast  the  boundary  value  problem  thus  far  established 
in  dimensionless  form,  and  for  this  purpose  introduce  the  following  dimen- 
sionless  parameters: 

x  =  x/a  ,      y=  y/a.  ,      z=z/a  }   d=  3/a    ,     k=k  /<z     ■> 


H=H/za}     h=k/a}      y^G-^oi/q     ,    t  =  crt 

<p  =     CT<J)/gd 


(7) 


where  o  =  2tt/T,  T  being  the  wave  period.     The  dimensionless  time  is  denoted 
by  t,  and  h     denotes  a  dimensionless  Bernoulli   constant.     As  indicated  in 
Equation   (7),  all   the  coordinates  and  length  scales  are  made  dimensionless 
with  the  characteristic  dimension  of  the  fixed  object,  a". 

Using  these  parameters,  Eqs.    (2-6)  which  define  the  boundary  value 
problem  may  be  rewritten  concisely  in  dimensionless  form  as: 


V^cx.lj, z,t)-o 


4>M<*,-M/*>  =  ° 


in  the  fluid 


(8) 
(9) 


*nCx*y/£/£^  =  °    on    SCX,4,Z;  =  0  (10) 

^,},zt)Kx^£)  +  4><-X'liL>i)}Cx,i.t)-<ptxl1lz,-t:)  +V}(x,z,t)   =o        (11) 

4g*i1,*i*>+^[<l>CKl,Z/±f+4><*b*A)Z+  <fyx.},*,*)]  +?(*,*,*)  =   h*  (12) 

Equations  (11)  and  (12)  are  applicable  outside  the  body  on  the  free  surface 
when  the  body  is  surface  piercing. 

PERTURBATION  PROCEDURE 

According  to  the  methods  of  perturbation  analysis,  the  potential 
function,  free  surface  elevation,  and  Bernoulli  constant  may  be  expanded, 
respectively,  in  the  small  parameter,  e,  as 

<p(x,y,z,t)=  €  4>cx,y,-zt£)    +-<£*Q(x,y,z,-£)  +  0<&)  (13) 

j(x,Z,-L)     =  €?((/,z^)  *  £z  JCx,*,*)    +    0(eV  (14) 


where  e  relates  to  the  wave  amplitude  in  a  yet  undetermined  manner. 

It  may  be  noted  that  in  Eqs.  (11)  and  (12)  the  surface  elevation,"^  , 
and  consequently,  the  small  parameter,  e,  does  not  appear  explicitly. 
Accordingly,  it  is  necessary  to  expand  the  particular  functions  involved 
so  as  to  form  functions  which  do  contain  e  explicitly.  For  example, 
substituting  Eq.  (14)  into  (13)  and  expanding  each  term  in  a  Taylor 


(15) 


series  about  e  =  0  gives: 
<pC^hZ/t)  =  e  <%(x,o,z,l)  +^[<%(x,o,Z,t)+?/*,2.d)  4>  (x,o,Z,*)]  +  0(6')  (16) 

Similar  results  can  be  obtained  for  the  derivatives  of  <j>  which  appear  in 
Eqs.    (11)  and   (12). 

Substituting  Eqs.    (13-16)  as  well   as  expressions  similar  to  Eq.    (16) 
for  the  derivatives  into  the  boundary-value  problem  given  in  Eqs.    (8-12), 
the  separate  boundary- value  problems  for  $lt  <j>2,  etc.   are  obtained.     The 
first-  and  second-order  problems  may  be  isolated  and  are  given,  respectively, 
by: 

First-order: 

V^x.y,?,*;  =o  (17) 


c£(x,-h, 2,^)^0  (18) 

i) 

(fKx.y.E,*)    =  O         on         Scx,lj;z)   =0  (19) 


<jjcx,o,2,o  -  y%j:x,ztt)   =  o  (20) 


3<x,z,0    ■*-   d>cx,o,z,0  -O  (21) 

1  't 

Second-order: 
\7*<fe(x,y,i,0  =r  O  (22) 


cj>y(x,-h,£,^)    =  O  (23) 


10 


^(X.L),"?,*)   =    O  on  S(x,y,i)=0 


(24) 


^Cx^.z.t)  -  V  ?z(x,2.,t)=  -  icx.z.-t)  tiix.o  ,z.,t)  +  tytix.o.i.^l^it)      (25) 

^Cx,z,-L)+.  4>(x,0,Z,£)=- 7^,2,6;  $0cto ,1,-L)  - -L [fyx.o.Z.if  (26) 

+  <j>(x,ot3L,t)     -h    <$>Cx,o,z,±)  J  +    h\ 


Both  <!>!  and  <j>2  are,   in  addition,  subject  to  a  suitable  radiation  condition 
which  limits  the  scattering  disturbance  to  outgoing  waves. 

A  great  deal   of  similarity  is  evident  between  the  first-  and  second- 
order  problems.     In  fact,  the  only  basic  difference  lies  in  the  right 
hand  side  of  the  free  surface  boundary  conditions;  in  contrast  to  the 
first-order  free  surface  boundary  conditions,  the  second-order  boundary 
conditions  are  not  homogeneous,  the  right  hand  sides  of  which  being 
dependent  on  the  first-order  potential. 

Before  proceeding  to  the  solution  of  the  boundary- value  problems 
developed,   advantage  may  be  taken  of  the  fact  that  the  first-order  poten- 
tial will  be  represented  by  a  periodic  function.     Accordingly,  we  may 
define  the  complex  potential   ui(x,y,z)  as 


4>cx,ytz,£)   =  ab   F&[i   U, (*,</,*;  el    J 


(27) 


where  Re  denotes  the  real  part,  and  a  =  2-rra/L,  a  being  the  characteristic 
dimension  of  the  body  and  L  the  wave  length.  The  symbol  b  denotes  an 
unknown  real  constant.  It  is  also  appropriate  and  common  practice  in 
linear  interaction  problems  to  express  the  potential  <$>i   as  the  sum, 


u 


4>   =  <P*    +   <jf  (28) 

where  $     denotes  the  incident  wave  potential,  and  $     denotes  the  scattering 
potential  which  is  due  to  the  presence  of  the  rigid  body.     The  complex 
potentials  u1.   and  u1  are  then  defined  in  relation  to  ^and  <j>&   according  to 
the  form  of  Eq.    (27)  so  that  it  follows  that 

u<    =    "*    +    U?  (29) 

It  is  recognized,  moreover,  that  the  incident  wave  potential,  $\t  must 
satisfy  the  first-order  problem  when  no  body  is  present,  i.e.,  <|>*  must 
satisfy  Eqs.    (17),   (18),   (20),  and  (21).     These  equations  represent  simply 
the  boundary-value  problem  for  the  well-known,   first-order  progressive 
wave.     The  solution  to  this  in  terms  of  the  present  notation  is  given  by: 

uf  =  --JL  cosMa(h+y)l     e/ax  (30a) 

In  Eq.  (30a)  a  is  defined  in  terms  of  v  by  the  familiar  expression  from 
linear  wave  theory 

V-  cr^a/g  =  a    fankiak)  (30b) 

where  h  denotes  the  dimensionless  water  depth. 

Now,  substituting  Eqs.    (30a),   (29),  and  (27)  into  the  first-order 
problem  defined  by  Eqs.    (17-21),  and  eliminating  n     between  Eqs.    (20)  and 
(21),  the  following  boundary-value  problem  is  established  for  the  first- 
order  complex  scattering  potential: 
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VzU^Cx,H,z)   =  o  (31) 


U*y(x,-M)  =  o  (32) 

U*C*.4,2)  = J i —   ns  sinK[o-Ch+9)]  +■  L  ftx  coshLaCK+y)]    eaX  (33) 

o/z      5c*,  y,z.)  =  O 

U?ifKo/z)-yU,sCx,o/z)   =0  (34) 


The  components  of  the  unit  vector  directed  normal   to  the  immersed 
surface  into  the  fluid  are  defined  as  n^n^j^fft^  .  In 

addition  to  Eqs.    (31-34),   u*(x,y,z)  is  also  subject  to  the  usual 
radiation  condition  which  allows  only  outgoing  scattered  waves  at  a 
great  distance  from  the  body. 

Before  dealing  with  the  solution  to  the  first-order  scattering  problem, 
the  second-order  problem  will   first  be  established.     Proceeding  in  this 
direction  it  is  appropriate  in  view  of  the  linearity  to  express  the  solu- 
tion to  the  second-order  problem,  Eqs.    (22-26),  as  the  sum 

*z   =    4  +    €  (35) 

where,  similar  to  the  first  order  problem,   <f2  denotes  the  second-order 

incident  wave  potential   and  4?  denotes  the  scattering  potential. 

We  again  envoke  the  condition  that  when  no  body  is  present,  there  is 

no  scattered  wave  so  that  <t>s  =  41s  =  0.     Moreover,  under  these  conditions 

1         2 

the  boundary  conditions  on  the  immersed  surface,  Eqs.    (19)  and  (24),  are 
not  applicable.     Accordingly,  we  may  substitute  <j>     =  4>1  and  <}>2  =  <j>   ,  along 
with  the  Eq.    (30)  and  (27),  into  Eqs.    (22),    (23),    (25),   and   (26)  and 
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obtain,   upon  eliminating  n  between  Eqs.    (25)  and   (26),  the  following 

boundary-value  problem  for  the  second-order  incident  wave  potential: 

Vz<$(x,Lilzl±>)    =     O  (36) 

<t}cx,-h,z,t)  =  o  (37) 


4>(x|0/2^;  +   y  <p\x,o,z,4)  --%  b2(al-  vL)  sm[z(ax-£ )]  (33) 

The  periodic  solution  to  the  boundary-value  problem  specified  by 
Eqs.    (36-38)  may  be  expressed  as 


cp2=  3  abz  y>  pc[i  Utz(xiH)  eUi] 

where  the  second-order  complex  potential,  u1 ,  is  given  by 


(39) 


ufoy)  =  -  T7T   coshUaOi+y)'}    <~zclK  (40) 

This  expression  is  familiar  in  wave  theory  and  is  of  exactly  the  same 
form  as  that  given  by  Eagleson  and  Dean  (1966)  for  second-order  Stokes' 
waves. 

To  complete  the  establishment  of  the  boundary-value  problem  for  the 
second-order  scattering  potential,  it  is  necessary  to  substitute  the 
forms  4.  =  t})1  +  <t>s   and  4>  =  <j>r  +  <$>s   into  Eqs.  (22-26)  and  eliminate  n 
between  Eqs.  (25)  and  (26).  In  addition,  the  known  expressions  for  the 
incident  wave  potentials,  Eqs.  (30)  with  Eqs.  (29)  and  (27)  for  the  first- 
order  potentials,  and  for  the  second-order  potentials,  Eqs.  (40)  with 
Eqs.  (39)  and  (35),  are  utilized  to  obtain: 
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V2<J>s<x,q,z,£  ;  =o  (41) 


<jf  Of.-M,*)   =■  o  (42) 

4>S(x,  cj  ,£,£)  =  -  4.  ^^_>i^ }{nxcoshL2ach+y)J-  Lrcy  s.nhL2aCh+y)]feLZCaX"t)  (43) 


(44) 


'£  'y 


In  Eq.    (44),   U(x,z)ns  a  time  independent  function  which  is  generated  by 
substitution  of  <j>     into  Eqs.    (25)  and  (26).     However,  since  the  function 
is  not  needed  in  the  second-order  theory  for  purposes  of  evaluating  the 
pressures  and  forces,   it  is  for  brevity's  sake  not  written  out  in  full 
here. 

The  boundary-value  problem  now  established  for  <j>     as  given  in  Eqs. 
(41-44)  is  linear,   and  except  for  U(x,z)  occurring  in  Eq.    (44),   is  time 

-l*2t 

dependent  like  e    .  Accordingly,  it  is  appropriate  to  express  the 
solution  in  the  form 

where  the  second  term  denotes  the  part  of  the  complex  potential  which 
is  independent  of  time.  By  substitution  of  Eq.  (45)  into  the  boundary- 
value  problem  described  by  Eq.  (41-44),  separate  boundary-value  problems 

s     s  s 

arise  for  u2  and  u2  •  However,  as  will  be  evident  subsequently,  $ 
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only  is  needed  in  order  to  evaluate  the  hydrodynamic  pressures  and 
resulting  forces  to  the  second  order  and,  consequently,  the  time 

c 

independent  part  of  <p     is  of  no  interest.  We,  therefore,  dismiss  it 
from  further  consideration  and  concentrate  on  the  solution  u2. 

Substituting  Eq.  (45)  into  Eq.  (41-44)  the  following  boundary 
value  problem  is  established  for  the  second-order  scattering  potential 


V^O.h.h)  =  o 


(46) 


Uly(x,-h,z)  =  O 


n  Smh    (nh)  l     J  -J 


zax 


r~ 


^  (X,o,z)  -    4-y  uf(x,o,z)  -    f(x,z) 


(47) 


(48) 


(49) 


where 


f  <*x,H)  =  |  A/i  «*   u;  *  «,'(W,;3  -  5  V  «,J  )  *-$.  «J  (  U^-VUfJ     (50) 


1  J  y=  o 


The  radiation  condition  limits  the  scattered  waves  to  outgoing 
regular  waves  and  is  given  by 


\im     V/c" '  (1 :;-   L  Qz)  US      =■ 


o 


(51) 


where  a2  denotes  the  dimensionless  wave  number  for  the  second-order 
problem  which  is  defined  by 


4-V  =  QL    -f-anh(Qji) 


(52) 
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INCIDENT  WAVE 

At  this  juncture  it  is  appropriate  to  completely  define  the  incident 

* 
wave  and  specify  the  unknown  constants  b  and  h2.  Solving  Eqs.  (21)  and 

(26)   for  t\     and  n   ,  respectively,  and  substituting  the  results  into  the 

expression  for  the  elevation  of  the  free  surface  defined  by  Eq.    (14) 

gives 

(53) 
ZvL    'x  y  Z  J  J 

Equation  (53)  expresses  the  free  surface  elevation  in  terms  of  the  total 
potentials  and,  consequently,  includes  the  effect  of  the  scattered  wave 
as  well  as  the  incident  wave.  However,  it  is  presently  of  interest  to 
obtain  an  expression  for  the  free  surface  elevation  of  the  incident  wave 
alone  in  the  absence  of  the  object.  For  this  purpose,  therefore,  we  set 

c  c 

$     =  $     =  0  and  evaluate  Eq.    (53)  using  the  known  expressions  for  the 
i  2 

incident  wave  potentials,   Eqs.    (30)  with   (27)  and  Eq.    (40)  with   (39). 
The  result  is: 


+  -a~ 


I      4-r 
sirxhVah)  •> 


(54) 


/'/ 


where  nr  denotes  the  dimensionless  surface  elevation  due  to  the  inci- 
dent wave  in  the  absence  of  the  body.  However,  the  constant  term  in 
Eq.  (54 )  must  vanish  so 

K    -  *&*  (55  ) 

The  remaining  terms  in  Eq.  (54)  represent  sinusoidal  variations,  the 
second-order  contribution  having  twice  the  frequency  of  the  first. 

If  we  define  the  wave  height  as  the  elevation  difference  between 
the  trough  and  crest  of  the  incident  wave,  then  eb  must  represent  the 
wave  amplitude  since  the  second-order  contribution  to  the  surface 
elevation  in  Eq.  (54  )  is  the  same  at  both  the  crest  and  trough,  the 
difference  being  zero.  Therefore,  in  terms  of  the  dimensionless  wave 
height  as  defined  in  Eq.  (7)  we  have: 

€  b  =    H  =    H  /Za  (56) 

where  H  denotes  the  elevation  difference  between  the  crest  and  trough. 
The  incident  wave  profile  is  then  specified  by  the  dimensionless 
surface  elevation  as: 

}\x,*,i)**H  cos(ax-i)  +  H^-%-    CoshCah)  sz+  coshCzah))  (57) 

*  sink  Cah) 

s    cosUcax-i.)]  +  OCH3) 

Equation  (57)  agrees  with  the  expression  for  the  second-order  Stokes' 
wave  given  by,  for  example,  Eagleson  and  Dean  (1966). 


PRESSURE,  FORCES  AND  MOMENT 

The  forces  and  moments  acting  on  the  immersed  surface  are  deter- 
mined by  carrying  out  surface  integrals  of  the  pressure.  For  this 
purpose  the  pressure  is  obtained  from  Bernoulli's  equation  as 

P<x,y,M;  =  -  f>  4>_  -  ±p[4>_z+-  $*  +  4>^]-pS9  +  P3h  *       (58) 

This  expression  may  be  cast  in  dimensionless  form  by  use  of  Eq.  (7), 
(13),  (14),  (15),  (27),  (35),  (45),  and  (56),  and  the  result  carried 
to  the  second-order  in  wave  height.  The  resulting  expression  is  then 

fU,^,i)  =  -n-fta   Re[u,  eu]-  H^faftQg-u^  «,«  _  a-       {59) 

u*)eui]+\alx\\  |u,/  +  [u.4|l  +  j^.,1 

LA.  J 

The  dimensionless  pressure  coefficient  expressed  by  Eq.  (59)  is 
defined  as 

-f   -  .  P_  (60) 

'         P3Q- 

in  which  P  denotes  the  fluid  pressure. 

The  first  term  in  Eq.  (59)  denotes  the  hydrostatic  pressure  as 
y  represents  the  dimensionless  depth  beneath  the  mean  free  surface. 
The  second  and  third  terms  are  harmonic,  the  second-order  part  having 
twice  the  frequency  of  the  first-order  part.  The  final  second -order 
term  in  Eq.  (59)  is  independent  of  time  and  gives  rise  to  time-average 
or  steady  state  components  of  force  and  moments. 

The  components  of  the  dimensionless  force  and  the  moment  vectors 
may  now  be  expressed  in  terms  of  integrals  of  the  pressure  over  the 
wetted  surface  area.  That  is,  we  may  write 
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CL(t)  =  ff  f  9.  ds     ,       i  -  ',*/•■■  ^  (61) 

s 

where  the  dimensionless  force  coefficients  associated  with  the 
x,  y,  and  2  directions  are  defined,  respectively,  as 


C/0=    U^  (62a) 

P9a 


PS  a* 


(62b) 


C3(i)=    -£J±>-  (62c) 


in  which  F x(t),  Fy(t),  and  Fz(t)  denote  the  three  components  of  force. 
Similarly,  the  moment  coefficients  are  defined  as 


P9<> 


x*  (62e) 


Cift)-  fga*  (6/f) 


in  which  Mx(t),  My(t),  and  Mz(t)  denote  the  three  components  of  the 
wave-induced  moment. 
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The  functions  g-j   are  defined  as 
%  =  n*    >  9t=   nH     >  33  *    n*  (63a) 

g4=  CcUy)  n£  -  zr?y   ?     35^nx-xni;   c^=  xny-<iy+d)nK  (63b) 

The  forces  and  moments  in  Eq.    (6  2)  are  defined  relative  to  a  coordinate 
system  parallel   to  0(x,y,z)  but  shifted  downward  a  distance  d  as  is 
evident  from  Eq.    (63b).     The  unit  vector  directed  normal   to  the  wetted 
surface  is  defined  by   n=TnK  +fn^    +  k  nz        and  dS  denotes  an 
element  of  surface  area  made  dimensionless  with  the  characteristic 
dimension  of  the  body  squared,  a2. 

We  may  make  an  initial  ordering  of  the  force  or  moment  coefficient 
as  given  in  Eq.  (61)  in  terms  of  the  small  parameter,  e,  or  equivalent 
H,  as 

c>M*-fftS  Si  dS  ~  H\\r,   3ids-"lfft>  WS-Oe**)      (64) 

where  the  pressure  coefficient  is  expressed  as 

^p  =  ^y   +    H    ^     +     A/^  +.    £T(/;*)  (65) 

and  the  coefficients,  p0,  Pi,  and  p2,  are  defined  in  an  obvious  manner 
by  comparison  of  Eq.  (65)  with  Eq.  (59). 

It  may  be  noted  at  this  point  that,  if  the  body  is  completely 
submerged,  the  surface  integrals  in  Eq.  (64)  refer  simply  to  the  actual 
surface  area  of  the  body  which  is,  of  course,  specified.  However,  if 
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the  object  is  surface  piercing,  S  denotes  the  actual  wetted  area  which 
is  dependent  on  the  instantaneous  surface  elevation  of  the  water  on  the 
body.  Accordingly,  it  is  necessary  to  express  the  limits  of  integration 

on  the  wetted  surface  S  in  terms  of  n  and  carry  out  the  expansion  retain- 

2 
ing  terms  with  coefficients  up  through  H  . 

For  this  purpose  the  dimension! ess  differential  surface  area  may  be 

written  as  dS  =  dc  dl  where  dc  denotes  a 

dimensionless  differential  arc  length  on  the  wetted  surface  in  the 

horizontal  plane,  and  dl  denotes  a  dimensionless  differential  length 

in  an  orthogonal  direction.  Thus,.   it  is  possible  to 

express  dl  in  terms  of  y  so  that 


d5  =    dc  di  =  ^d^Ju. 


Ml 


(56) 


n, 


The  integrals  indicated  in  Eq.  (64)  are  now  carried  out  over  the  wetted 
surface  with  y  running  from  y  =  -e  to  y  =  n.  Accordingly,  a  typical 
integral  in  Eq.  (64)  may  be  expressed  in  the  form 

s  * 

y=  -e 

where  e  denotes  the  lower  limit  of  the  immersed  surface. 
The  integrals  of  the  type  given  in  Eq.  (67)  are  next  expanded  in  a 
Taylor  series  about  H  =  0  using  Leibniz's  formula  for  purposes  of 
evaluating  the  derivatives  of  the  integrals.  Carrying  out  these 

expansions,  substituting  the  results  into  Eq.  (64),  and  retaining  terms 

2 

up  to  order  H  yields: 


22 


CL(i)  =  -f[  H  9LdS  ~  "JJ  f<   ft  dS  '  "7//^  9:  ds 


(68) 


•f- 


where  S0  denotes  the  surface  area  below  the  plane  y  =  0  and  C0  denotes 
the  closed  waterline  curve  formed  by  the  intersection  of  the  y  =  0 
plane  with  the  surface  of  the  rigid  object. 

The  first  term  in  Eq.    (68)  represents  the  hydrostatic  force  or 
moment  on  the  body  where  the  associated  displaced  volume  is  defined  as 
that  beneath  the  y  =  0  plane.     However,  since  our  interest  is  in  the 
dynamic  forces  and  moments,  we  may  disregard  this  first  term  with  the 
understanding  that  the  buoyant  force  defined  as  such  should  be  accounted 
for  in  order  to  determine  the  final   total   force  or  moment. 

Finally,  applying  Eq.    (21}  and  (27)  as  well   as  the  definitions 
for  p0,  Pi  and  p2  we  obtain: 

Ci(t)±H-af[ReLu,eri*lgidS+Ht(:!gffRt[(*gul-U.) 


o 


-^//[KI^W,^  K|%  ^:-i]g;ds) 


4-      O(H') 
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Equation  (69>a)  may  also  be  written  in  coefficient  form  as 

$(*>  =  H  /T.   ^>sC8,l  -  t)+Hz(£.cosC£z-z+)+  tt**J  +  0(H3)  (69b) 

where  the  first-  and  second-order  force  (or  moment)  coefficients  and 
phase  shift  angles  are  defined  by  comparison  of  Eq.  §9a)  with  $9  b) 
as  follows: 

f[  el  ''  =   a/fa,  g.  c/s  (69c) 

z  RzYT-  n 


=r    -  /  )   U,     Q.    dC 


where  the  dimensionless  force  coefficients,   Fi.   and  F2.,  are  real. 

i     m 

The  steady  state  (nonperiodic)  force  (or  moment)  coefficient  is  also  defined 
by  comparison  of  Eqs.  (69  a)  and  (69b)  as 


*"/*-  £l(*rkw  -  0  ^  ^  dc 


(69  e) 


£.ff[ja^lu,^ju^  +  ±:-i]^6S 


-h 

4- 


The  first  term  in  Eq.  (69  b)  represents  the  linear  solution  having 
the  fundamental  frequency,  o.  The  next  term  which  represents  a  second- 
order  contribution  to  the  force  or  moment  is  also  harmonic  but  at  twice 
the  fundamental  frequency.  The  last  second-order  term  represents  the 
steady  state  contribution  which  is  independent  of  time.  This  latter 


24 


force  is  generally  referred  to  as  a  drift  force  in  ship  hydrodynamics. 
(See,  for  example,  Maruo  (I960).)  However,  drift  force  is  generally 
calculated  by  applying  the  momentum  equation  to  the  diffracted  wave 
rather  than  by  using  near  field  results. 

The  evaluation  of  the  force  or  moment  coefficients  defined  by  Eq. 
(39c-e)  represents  the  primary  object  of  this  paper.  However,  in 
order  to  carry  out  the  evaluations  indicated,  it  is  clear  that,  given 
the  geometry  of  the  rigid  immersed  object,  it  is  necessary  to  evaluate 
uL   and  u2  as  well  as  the  derivatives  of  ui  on  the  surface  area  denoted 
by  S  and  along  the  water! ine  curve  C  . 

SOLUTION  TO  THE  FIRST-ORDER  PROBLEM 

Having  now  established  the  need  for  the  potentials,  ut  and  u2,  we 
return  to  the  consideration  of  the  solution  to  the  boundary-value  prob- 
lem developed.     The  practical  method  of  solution  of  the  first-order 
problem  given  in  Eqs.    (31-34)   is  discussed  in  depth  by  Garrison(1974,1978) 
and,  therefore,  only  the  major  steps  will   be  outlined  here.     Following 
the  Green's  function  method  of  solution  we  write  uf,  as  the  integral 
over  the  immersed  surface  S     as: 

l/;(x,ytZ)  =-L  [($(5,1,1)    G<*.H,*'jZ,hSj-»)    dS  (70) 

n~TT  •'-'c 

where  U,n,e)  denotes  points  on  the  immersed  surface,  fx  U,n,c)  denotes 
the  unknown  source  strength  and  dS  =  dS/a2  denotes  the  differential 
surface  area  on  the  immersed  surface  made  dimensionless  with  the  charac- 
teristic dimension  of  the  object,  a.  The  function,  G  (x,y,z;£,n,C»v), 
denotes  the  Green's  function  which  must  satisfy  the  equation 
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VZG(x,yi*;Itl,S:  V)   =  &(*-$)  S(y--j)  S<«-l)  (71) 

in  which  6  denotes  the  Derac  delta  function.  The  function,  G(x,y,z;s,n,s*,v) 
must  also  satisfy  Eqs.  (32),  (34),  and  the  radiation  condition.  Such  a 
function  is  given  by  Wehausen  and  Laitone  (1960)  as 


f^K  /&* ^ e  f* ^a r** yl co sh &<*+*)]  j-         ,  (72a) 

7   -*  s,nh(Mk)   -  v  cosh(«h)  M 


in  which 


(72b) 


(72c) 


^  ^Ifr-If+Cz-D*]*  (72d) 


The  symbol,  aj,  is  defined  in  terms  of  h  and  v  as  the  solution 
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of  the  equation 

a,  -/-an/iCqh)   -  v  =  o  (72G) 

In  view  of  Eqs.    (72e)       and  (30b)   the  symbol   a  is  clearly  equivalent 
to  a.     However,   in  the  case  of  the  second-order  Green's  function,  the 
equivalence  does  not  hold  and,  therefore,  separate  notation  is  maintained, 
In  Eq.    (72a),  P.   V.   denotes  principal   value  of  the  integral. 

An  alternate  series  form  of  the  Green's  function  is  also  given  by 
Wehausen  and  Laitone   (1960)  as 

(73a) 


^(^fi+yh-v) 


k-~ 


where  J  and  YQ  denote,  respectively,  Bessel  functions  of  the  first  and 
second  kind  of  order  zero,  and  K0  denotes  the  modified  Bessel  function 
of  the  second  kind  of  order  zero.  In  Eq.  (73a)  aj  is  defined  by  Eq. 
(72e)  and  the  Quantities^  are  defined  as  real  positive  roots  of  the 
equation 


^k    ran  (^h)    -/-  V  =■  q 


(73b) 


The  solution  to  the  boundary-value  problem  stated  in  Eqs. (31 -34) 
is  given  by  Eq.(7^);  it  remains,  however,  to  determine  the  source  strength 
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function  in  order  to  evaluate  the  potential.  This  is  accomplished  by 
application  of  the  boundary  condition,  Eq.  (33),  which  results  in  the 
following  integral  equation  from  which  fx  may  be  determined: 

^/fa/.X>|S<*.S.^^^  (74) 

So 

The  normal  derivative  of  G  evaluated  on  the  immersed  surface,  as  required 
in  Eq.  (74),  is  determined  in  a  straightforward  manner  by  differentiation 
of  either  Eq.    (72a)  or  (73a). 

SOLUTION  TO  SECOND  ORDER  PROBLEM 

The  solution  to  the  second-order  problem  may  be  formulated  by  application 
of  Green's  theorem  to  the  fluid  region  denoted  by  -J^in  Figure  2.     Applying 
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FIGURE  2   Region  for  application  of  Green's  theorem. 


Green's  theorem  with  subjects  u^  and  G2  (the  second-order  Green's  function) 


and  assuming  that  G~  satisfies  the  equation, 


VzGt(XtytZ}St1,S)  -    &Cx~s>&(y-})Sce-s) 


(75) 


as  well  as  the  kinematic  condition  on  S.  given  by  Eq. (47) ,  the  radiation 
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condition  given  by  Eq.(51)  on    S^    and  the  homogeneous  free  surface  con- 
dition, 

G*/K*,*;S,TtS>  -4»Gz(*AZ;r,j,S)  -  &  (76) 

gives  the  result: 

«*&#*>=£  flatter)  <*„K%*;fJ,TJ-  G(*ti*;UJ)K<W<frs 

S°  (77) 

The  first  term  in  Eq . (77)  may  be  considered  to  represent  a  source  and 
doublet  distribution  over  the  immersed  surface  of  the  body.  However, 
Lamb  (1932)  has  shown  for  the  case  of  an  infinite  fluid  that  the  potential 
may  be  represented  as  either  a  source  distribution  only,  a  doublet  distri- 
bution only,  or  some  combination  of  both.  This  proof  is  easily  extended  to 
the  present  case  allowing  the  first  term  to  be  written  as  a  distribution  of 
sources  only.  Then,  if  the  free  surface  boundary  condition,  Eq.(49),  is 
used  in  the  second  integral  the  result  is 

(78) 

-&$ Win  <**<*«•*'' w>  ds 


Ss 


in  which  f±(f,?fT)       is  the  unknown  source  strength  distribution. 

The  second-order  scatter  potential  defined  by  Eq. (78)  satisfies  Eqs. 
(46,  47,  49-51).  The  remaining  kinematic  boundary  condition  on  the  surface 
of  the  body  results  in  the  integral  equation 
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(79) 


Smh+(4h) 


+  ^ff  fa?)  ^(*<y<*;!>°<*;+y<«*>  ^ 


'Ss 


which  is  to  be  satisfied  for  x,y,z  on  S0    . 

In  both  Eq. (73)  and  (79)  a  practical  interpretation  of  the  limits  of 
the  free  surface  area,  Ss  >  should  be  understood.  In  deriving  Eq.(78) 
through  application  of  the  Green's  theorem  the  outer  limit  of  Ss     was 
considered  to  be  large  enough  that  the  radiation  condition  was  satisfied 
on   S,*,   .  in  practical  application  this  occurs  rather  rapidly. 

In  view  of  the  conditions  which  were  placed  on  the  Green's  function  in 
the  above  development  G~  is  clearly  \/ery   similar  to  the  Green's  function 
used  in  formulating  the  first-order  problem.  In  fact,  if  V  is  replaced 
by  4V  and   #.   is  replaced  by  Ct,_    in  either  Eq .  (72)  or  (73)  the  re- 
sulting function  is  the  appropriate  Green's  function  for  solution  of  the 
second-order  problem.  That  is, 

Git**,*;*,!,-*)   =  Gtov,*;  5,1,3;+*,**)  (so) 

where  G  is   given  by  Eq . (72)   and   (73)   and  0Lis  defined  through  Eq.(72e) 

az    -Unk  azh     -    4V    =  O  (81a) 

In  the  case  of  the  alternate  series  form  of  G2  given  in  Eq. ( 73) ,  the  roots 
are  defined  through  Eq. ( 73b)  as 


^k     ^n^h)      +     +V    =   O 


(81b) 
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An  alternate  form  of  Eq.(72a)  for  the  sources  on  free  surfaces  is  given  by 


4-yh  +  s/n/?2-^/? 


(82) 


This  if  the  form  of  the  Green's  function  which  is  obtained  by  setting  |so 
in  Eq.(72a)  and  which  was  also  given  by  Wehausen  and  Laitone  (1960)  for 
the  case  of  a  harmonic  pressure  distribution  on  the  free  surface.  The 
nonhomogeneous  form  of  the  free  surface  boundary  condition  given  in  Eq.(49) 
is  equivalent  to  the  boundary  condition  for  the  first-order  problem  for  a 
pressure  distribution  of  the  free  surface. 

NUMERICAL  PROCEDURE 

Eq.(69)  clearly  indicates  that  the  determination  of  the  forces  and 
moment  acting  on  the  immersed  object  rests  on  the  determination  of  the 
potentials  u,-,  u~  and  derivatives  of  u-,  at  points  on  the  immersed  surface. 
The  first-order  scattering  potential  is  specified  by  Eq.(70)  in  terms  of  the 
first-order  source  strength  function,  f , .  Thus,  to  determine  u-, ,  it  is 
first  necessary  to  solve  the  integral  equation,  Eq.(74),  for  f -, . 

We  may  develop  a  numerical  solution  to  Eq . (74)  beginning  with  the 
partitioning  of  the  immersed  surface  S  into  N  subdivisions  or  facets  of 
area  AS;  ,  each  with  a  nodal  point  at  its  center  located  at  the  point 
(x.,  v.,  z.).  Recognizing,  moreover,  that  f,  (|,?,J)  is  a  well-behaved 

vj       U       J 

function  for  smooth  bodies  we  may  define 

k  =  ssj7sji/"iv*'*/iy s'ntt«(Hi)kin>iw&C6*fa(i»'/^  (84) 
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fij    =  $(*,%,*;> 


(85) 


and  accordingly  we  may  approximate  Eq.(74)  by  the  complex  matrix  equation: 

(XLj(y)  £■    -    hL  i,J*i,at'--N  (8f>) 

Once  Q{.  .(^)    has  been  evaluated  by  use  of  Eq.(83)  the  solution  of  Eq.(86) 
may  be  carried  out  to  determine  f. .  at  each  nodal  point  on  the  surface  of 
the  body.  A  discussion  of  the  details  of  the  numerical  evaluation  of  C{.  ; 
including  the  evaluation  of  the  i/R  -  singularity,  has  been  given  by 
Garrison  (1973). 

Having  determined  the  source  strength  function,  the  first-order  poten- 
tial and  its  derivatives  on  the  surface  of  the  object  may  be  determined  from 
Eq.(70).    For  this  purpose  we  may  replace  the  surface  integrals  with 
summations,  writing 

of.  =  fi.t»)  f.  ,         l,j  =,iZi...  A/  (87) 

U'H  =  P*i}y)    tj       ,      i.  '*  '<*,•  ■  «  (33) 

Uf>r^y)f'J       '        *"'"'<</~"  (89) 

U\t  =  A4-5v;    $J   i  l'-'J~  '>*,*-"  (90) 


,      c  C  ■*-  U 

where  u1   ,   u:      ,  etc.   denote  functions  evaluated  at  the  i       nodal   point 
i  X, 

on  the  imiTiersed  surface,  S0.     The  complex  matrices  occurring  in  Eqs. 
(37  -90  )  are  given  by 


#/v;  -■-  ^rjf  Gte ,  »,£■ ;  I  h  T;  y)  ^S  (gi  } 

AS- 
j 

fii<y)    ■  :fk?[[ ^(X<^'^^U;^  te  (9?) 

as 


(93) 


P*oly)  '  <&fj[fi<*>  •">*>*'.  r;  y>  ^  {9A ) 


In  evaluation  of  the  integrals  in  Eqs.(83)  and  (91-98)  it  is  generally 
adequate  to  simply  evaluate  G  and  its  derivatives  at  the  centroid  of  the 
panel  and  multiply  by  its  area.  However,  when  point  i  is  either  equal  to 
j  or  near  j  it  is  necessary  to  take  more  care  in  integrating  the  1/R- 
singularity  in  G.  The  procedures  given  by  Garrison  (1978)  were  used  to 
evaluate  the  integrals  on  quadralateral  panels. 

To  obtain  a  solution  to  the  second-order  problem  it  is  necessary  to 
first  evaluate  the  function  f  (j,3)  required  in  Eq. (49)  and  as  defined  by 
Eq.(50).  For  this  purpose  we  proceed  again  numerically  by  dividing  the 
mean  free  surface  (or  y  =  0  plane)  in  the  vacinity  of  the  body  into  K  area 
elements  and  evaluate  f  at  the  nodal  points  of  these  elemental  areas.  The 
numerical  procedure  is  based  on  the  use  of  Eq.s(87-94)  with  y-  =  0  for 
purposes  of  determination  of  u,   and  its  derivatives  at  points  on  the  free 
surface.  The  function  f  is  then  evaluated  at  the  nodal  points  by  use  of 
Eq.(50). 

The  solution  for  the  second-order  source  strength  is  obtained  numerically 
through  the  integral  equation,  Eq. (79) .  That  is,  Eq. (79)  is  written  as 


ocLim  £j  =  K-L  +q(+v>  £ 


t  =  /,  Z,  •  •  •  A/ 


(95; 


where  CY^ y  (<f  v)    is  defined  by  replacing  V  by  ^-y  in  Eq.(S3)  and 

hz.  =  — J      .  //?y;  stoAb<itt*ft)]  1 1  I7K£  coshl£4(h*ft)]]el      L       ( 96 } 

Once  the  (X  -matrices  in  (9b)  hove  been  evaluated  it  may  be  solved  for  the 
source  strength  distribution  on  the  immersed  surface. 
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The  hydrodynamic  pressure  and  resulting  forces  carried  to  the  second 
order  as  given  in  Eq . (59)  and  (69),  respectively,  require  only  the  second- 
order  potential  rather  than  its  derivatives.  Once  f~  is  known,  u~  may  be 
obtained  through  Eq.(78)  which,  when  written  in  indicial  form,  becomes 


<  -  A/^  &  -  fik  (*»>  C    LZ'X'Z 


(97) 


where    fiijC+y)  is  defined  by  Eq.  (91 ) 


NUMERICAL  RESULTS 

A  computer  code  has  been  developed  for  the  special  case  of  a  vertical 
circular  cylinder  (pile).  MacCamy  and  Fuch's  solution  was  used  for  the 
first-order  solution  and  for  purposes  of  evaluating  f  .  The  second-order 
solution  was  calculated  on  the  basis  of  the  distributed  singularities  as 
given  by  Eq. (78). 

The  particular  example  calculation  was  carried  out  for  a  cylinder  placed 
in  water  one  radii  deep.  The  dimensionless  frequency  parameter  was  varied 
over  the  range  0<  V<   1.2.  The  results  for  the  first-order  horizontal  force 
coefficient,  the  second-order  horizontal  force  coefficient,  steady-state 
force  coefficient  and  phase  angles  is  shown  in  Figure  3.  The  results  appear 
to  indicate  that  the  second-order  effects  are  greatest  at  the  lower  freq- 
uency range,  i.e.,  for   v><0.6. 

In  Figure  4  and  5  the  dimensionless  horizontal  force  is  plotted  for  a 
complete  wave  cycle  for  two  different  frequencies  (wave  lengths).  These 
results  which  correspond  to  rather  steep  wave  show  a  rather  sizeable  non- 
linear effect.  Also,  the  effect  of  including  second-order  effects  is  to 
shift  the  maximum  force  towards  the  phase  of  the  wave  crest. 
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